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Summary 


This report discusses activities under NASA Contract NCR 39-009-077 for 
the period January 1, 1975 to June 30, 1975. During this period, the work was 
divided Into two phases; the results for each phase may be summarized as follows: 

1. Oscillatory Combustion of Liquid Monopropellant Droplets . Major emphasis 
was placed on this portion of the Investigation. Analytical efforts were 
continued for the prediction of open-loop droplet response. The zero-order 
calculations required In order to predict response were completed for 

a pressure range of 1-100 atm, over the range of droplet sizes Including 
small droplets which evaporate with little reaction effects and large 
droplets where the reaction zone Is very close to the surface of the 
droplet. Predictions of gasification rates and liquid surface temperatures 
were compared to earlier monopropellant droplet combustion measurements 
In a combustion gas environment. Good agreement was obtained between 
theory and experiment over the range of the measurements. 

Since the zero-order results were satisfactory, effort has turned to 
response calculations using the first order analysis. Calculations at the 
large droplet limit were completed with both quasisteady and unsteady 
. gas models for pressures of 1-100 atm. The open-loop droplet response 
factor was greater than unity for frequencies extending beyond the charact- 
eristic frequency for the liquid phase thermal wave. A second peak Is 
found near the characteristic frequency for gas phase transients, but 
response In this region cannot be Interpreted from simple response factor 
considerations. Response at the evaporative limit Is negative at low 
frequencies, peaking with a positive value near the liquid phase character- 
istic frequency. The maximum response factor at this limit Is less than 
unity. Calculations are In progress In the Intermediate regime between 
these limits and these results will be reported later. 

Consideration of droplet response at low frequencies, where the 
droplet size changes appreciably during a period of oscillation, will 
be undertaken in the next report period using the zero-order analysis 
for gasification rate predictions. The related problem of determining 
spray response from a knowledge of droplet response, at higher frequencies, 
will also be undertaken at that time. 

2. High Pressure Bipropellant Combustion . Work in this area involved extending 
earlier bipropellant combustion models, valid at low and high pressures 
respectively, by developing an intermediate pressure model which provides 

an Improved treatment of polar compounds. A second phase of this work in- 
volves a generalized analysis of high pressure combustion which allows a 
simple determination of the conditions required for critical combustion. 
Difficulties have been encountered with the current Intermediate pressure 
model as its high pressure limit is approached. Work is in progress to 
resolve this problem. 
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Semi-Annual Report for 1975 on the Investigation 
of Critical Burning of Fuel Droplets 


1. Introduction 


The objective of this investigation was to continue earlier work on the 
steady and unsteady combustion of liquid fuel droplets under rocket engine conditions. 
Emphasis has been placed on consideration of combustion at elevated pressures 
and temperatures, representative of typical combustion chamber environments. 

The results of this study have technical application to the design of liquid 
fueled rocket engines and the determination of the combustion instability 
characteristics of these engines. 

This report gives a summary of progress on the investigation for the period 
January 1, 1975, to June 30, 1975. During this report period, the work was divided 
into two phases, as follows: 

1. Combustion response of monopropellant droplets. This portion of the 
study is concerned with extending the analysis of the "open loop" response 
of burning monopropellant strands (1,2), to the case of droplets. The 
basic method employs a perturbation analysis allowing for both liquid 

and gas phase transient effects. Aside from the overall general problem, 
several asymptotic regimes can be delineated. This includes evaporation 
without reaction, a quasisteady gas phase, and the highly reactive regime 
where droplet curvature effects are negligible. A second category involves 
the low frequency limit where a single oscillation comprised a large 
fraction of the lifetime of a droplet and different procedures must be 
used to compute the response. 

2, High Pressure Droplet Burning in Flowing Combustion Gases. The objective 
of this phase of the investigation is to extend earlier work on bi- 
propellant droplet combustion (3-6) . An intermediate pressure model is 
under development to provide a continuous transition between present low 
and high pressure theoretical models; a generalized analysis is also being 
undertaken in order to obtain guidelines on the characteristics of liquid 
droplet combustion at elevated pressures; and finally, these results are 
being integrated in order to determine spray combustion characteristics at 
supercritical conditions for comparison with experimental results. 

During the present report period, major effort was devoted to the area of 
droplet response and activities in this area will be described in some 
detail. A brief description of work in the area of high pressure droplet 
combustion will also be given. 

Reports and papers issued during the period of this report are as follows: 

G. M. Faeth and S. Chanln, "1974 Annual Report on the Investigation of 
Critical Burning of Fuel Droplets," Contract NGR 39-009-077, The Pennsylvania 
State University, University Park, PA, February 1975. 
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2. Combustion Response of Monopropellant Droplets 
2.1 Introduction 


The objective of this phase of the investigation involves the determination 
of the open-loop combustion response of monopropellant droplets to imposed 
pressure oscillations. Knowledge of the frequency ranges where the combustion 
response is a maximum, allows the designer either to adjust geometry in order to 
avoid having characteristic chamber frequencies in this region or provide damping 
tuned to this frequency range. 

The present work is a continuation of earlier efforts on burning liquid mono- 
propellant strands (1,2). In the earlier work, a direct experimental confirmation 
of the response analysis was obtained. Using these results, the present extension 
to the case of droplets can be approached with some confidence. 

The theoretical approach involves the use of perturbation analysis. The 
Imposed pressure field is decomposed into a constant mean pressure and an os- 
cillatory component having a small amplitude, e. Decomposing all other quantities 
in the same manner, solution of the zero order ( e°) conservation equations give 
the steady state combustion characteristics. Solution of the first order (e) 
conservation equations provide the desired response characteristics; largely 
manifested by m*, the amplitude of the gasification rate perturbation (an exception 
to this occurs in the case of an unsteady gas phase as noted later) . 

Within the perturbation analysis, several asymptotic regimes may be delineated 
where substantial simplifications are available for the solution. These regimes 
are considered separately in the analysis in order to reduce computational effort 
while still covering a broad range of system variables. Solutions in the asymptotic 
regimes are also valuable for checking the general analyses which is valid over the 
entire range. 

The asymptotic regimes amy be subdivided with regard to reactivity (droplet 
size) and frequency, as follows: 

1. Reactivity. For a given mean pressure, the distance between the droplet 
and the zero order reaction zone varies with the size of the droplet. For 
very large droplets, the reaction zone is close to the droplet and 
curvature effects are negligible (7), allowing direct use of the earlier 
strand combustion results (1,2) for the estimation of response (in the 
following this regime is called the large droplet limit) . For very 
small droplets, the zero order reaction zone is located far from the 
droplet and fuel simply evaporates and diffuses away from the droplet 
with little reaction. In this region, the analysis may be simplified by 
neglecting reaction (small droplet limit) . Between these limits there is 
an intermediate size range where both reaction and curvature are important; 
here the complete solution must be considered. 

2. Frequency. At high frequencies, unsteady effects are important in both 
gas and liquid phases and the general analysis must be used (an exception 
occurs at very high frequencies where the liquid no longer responds and 
only the gas phase is unsteady, but this regime is not considered here). 

As the frequency of the pressure oscillation is reduced, a region is 
approached where the gas phase adjusts rapidly to transient effects, 
while the liquid phase is still unsteady due to its lower thermal diffusi- 
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vity. In this region, the gas phase is quasisteady, providing a sub- 
stantial reduction in the complexity of the calculations. 

At very low frequencies, the time of an oscillation is no longer 
small in comparison to the droplet lifetime. In this region the present 
analysis is no longer valid, and procedures more analogous to vaporization 
response calculations (8, 9, 10) are more appropriate. This regime 
will be considered after the calculations for the perturbation analysis 
are concluded. 


2.2 Basic Formulation 

A discussion of the assumptions of the analysis; the derivation of the perturbed 
conservation equations; and the notation of the analysis are all given in Ref. 7. 

In order to avoid duplication, these results will not be repeated here and the 
discussion will continue from Section 3.3.4 of Ref. 7. 

In the perturbation analysis, any generic dependent variable, G, is assumed 
to have the form 

G (r, t) = G^ (r) + e G^^ (r) e^^ (1) 

where e represents the small amplitude of the pressure oscillation and all quanti- 
ties are dimensionless. 

Substituting these quantities into the transient conservation equations and 
separating into like powers of e gives the zero and first order equations which 
must be solved. 


2.3 Zero-Order Solution 


2.3.1 Zero-Order Equations 

The zero order, steady state, equations are as follows (7); 

m, = m = const, 
fo o 


In the liquid phase, r < 1 

d 


, dT m 

, T ) = 0 

dr dr Pf°f ° 


In the gas phase, r > 1 


d0 


5-(r^T ^-m0) = O 

dr ' o dr o o' 

dT 

~ (r^ T ■— - m T ) + r^ q^"“ A T (0 - T exp (- E/T ) 

dr o dr o o o ' o o • o 


= 0 


( 2 ) 


(3) 


(4) 

(5) 
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The boundary conditions are; 
r = 0 
r = 1, 


'^o = ^i 


( 6 ) 


dT 


dT 


s- s+ 


(7) 


d (6 - T ) 


T„ (— (6^ - - q) 

3 + 


dr 


( 8 ) 


0o - To = a q exp (- L^/T^) 


r = oo 0=T^ = 1 

o o . 


(9) 

( 10 ) 


These equations are sufficiently general to allow for any ambient conditions, 
either adiabatic or non- adiabatic . Since most combustion chambers are essentially 
adiabatic, and monopropellants do not involve a local enthalpy defect in any stream, 
(This is only strictly true under the conditions of the present model which allows 
for temperature gradients in the liquid phase. If there is bulk heating of the 
liquid small enthalpy defects will arise in the spray.) a condition of adiabatic 
combustion can be applied. This leads to the zero order compatabllity condition 
between dimensionless properties as follows: 

q = 1 - 8T^ + L. (11) 

When Eq. (11) is placed in dimensional form, it provides a computation for the 
flame temperature. 


2.3.2 Zero-Order Liquid Phase 

The solution of Eq. (3) can be obtained immediately for the boundary condi- 
tions T = T, at r = 0 and T = T at r = 1, as follows: 

o 1 o os 


T - T, m 

o i f o 

T - T 

os i ^f f 


(1 - i)> 


( 12 ) 


Evaluating the derivative of Eq. (12) at r = 1 and substituting into Eq. (7) 
yields the boundary condition at r = 1, 

dT 

= *0 - B Tj + L) 

s+ 


(13) 
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2.3.3 General Zero-Order Gas Phase 

For the case of adiabatic combustion. Integration of Eq. (4) once, 
application of the boundary conditions at the liquid surface, and the 
adiabatic flame condition of Eq. (11), yields 

, de 

r T ^ - m 9 = - m 

o dr o o o 


( 14 ) 


Integrating Eq. (14) again, after transformation to eliminate explicit dependence 
on T^, yields after application of the outer boundary condition 

6„-l (15) 

as a consequence of adiabatic combustion. Substituting Eq. (15) Into Eq. (8) 
yields 


dT 

(T = m (q + T - L) 

odr'. o^ o 's 
s+ 


(16) 


= 1 - a q exp (- L^/T^^) 


(17) 


which provides a transcendental equation which may be solved for the temperature 
at the liquid surface. The outer boundary condition, Eq. (10) is unchanged. 

The zero order problem is now reduced to the solution of a single, nonlinear, 
second order, ordinary differential equation, Eq. (5). Equations (10), (16) 
and (17) provide three boundary conditions to define the problem and determine 
the eigenvalue, m^. 

In the general case, Eq. (5) is nonlinear and must be integrated numerically. 
Since any practical range of numerical integration must be finite with respect 
to r, the outer boundary condition, Eq. (10), was handled by considering the 
asymptotic behavior of Eq. (5) for large r. 


To find the form of Eq. (5) for large r, let 

T = 1 - C T 
o o 

where ^ is a small parameter. The radial coordinate is also stretched by 
defining 

j.a 

r = ^ n 


(18) 


(19) 


where a is a constant to be determined by proper matching, 
lowest order terms, the equation becomes: 


Keeping the 


d o . . ^ o 

dri dq ^ ®o ^ dq 


_2a +n- 12. . „ (1) n 

- i q A exp (-E) T^' 


= 0 


( 20 ) 
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The form of the asymptotic solution depends upon the reaction order, n. 

The combustion of hydrazine which was considered in Refs. 1 and 2, was best 
represented by n = 2, which is also generally true for other monopropellants. 

In this case, a proper match of terms can only be achieved by setting a = 0, 
and balancing convection against diffusion. The resulting asymptotic equation is: 

( 1 ) 


^ 7 dT 

d , 2 o 

d? dT 


( 1 ) 


- m T 
o o 


) = 0 


( 21 ) 


with the outer boundary condition T^=0,r=«>. Therefore, unlike the strand 
combustion case (1,2), the asymptotic solution for droplets for n - 2 does not 
Involve the reaction rate term. This occurs since the area for heat conduction 
increased as r^ rather than remaining constant, increasing the significance of heat 
conduction. 


The solution of Eq. (21) satisfying the outer boundary condition yields: 


T^ = 1- C [1 - exp (-m^/r)] (22) 

at large r, where the small parameter ^ has been absorbed into the unknown 
constant, C. The constant C must be determined by matching Eq. (22) with the 
numerical solution of Eq. (5). 

Given A, the solution procedure involves guessing a value for C and m^. 
Starting values for the numerical solution . determined from Eq. (22) 

^^o’ Hr '' 

at some large but finite value of r. Equation (5) is then integrated to the 
droplet surface, r = 1 using the fourth order Hamming predictor-corrector method. 

At r = 1, Eqs. (16) and (17) must be satisfied, determining and C uniquely. 

Since Eq. (5) is nonlinear, a double iteration technique must be used. The present 
procedure employed the Newton-Rapheson method for this iteration. Once m and C 
are determined, the outer boundary condition, Eq. (10), is only satisifed at the 
true mathematical infinity, through Eq. (22). 

Initial work on the solution of Eq. (5) involved a "shooting technique" 
Integrating out from the droplet surface and applying the boundary condition of 
Eq. (10) at a finite value of r. This procedure converged more slowly than the 
asymptotic method, and gave slightly higher values for m . The "shooting method" 
became virtually unworkable at large values of A where tRe reaction zone is very 
close to the droplet. This experience is generally similar to that encountered 
during the investigation of strands reported in Refs. 1 and 2. 


2.3.4 Gas Phase Solution for Small A. 


For a second order reaction, the parameter A in Eq. (5) is proportional to 
(p^ • Therefore, for small droplet sizes, A approaches zero and the reaction 

terms may be neglected in the conservation equations. Under these conditions, 

Eq. (5) may be Integrated to yield the following implicit equation for the 
temperature distribution. 


m 

o 





(23) 
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P 


The gasification rate, m , may be determined by applying the outer boundary 
condition, =» 1 when r = ", in Eq. (23) to yield 

“o = ^ ■ ^08 [1 - (1 - (24) 

The liquid surface temperature, can be determined from Eq. (17). 


2.3.5 Zero-Order Gas Phase for Large A 

At any pressure, A increases as the size of the droplet increases. With 
increasing A, the zero-order reaction zone moves closer to the droplet, reducing 
the Influence of curvature on the solution. This is also accompanied by an 
increase in m^, although T^^ remains unchanged. Reference to the liquid phase 
solution, Eq. (12), indicates a rapid decay of temperature as one moves into the 
liquid phase, reducing the effect of curvature in this region as well. 

The above discussion indicates that at large A, the process approaches a 
condition where temperature variation and reaction effects are confined to a 
thin region surrounding the liquid surface, with little effect of curvature in 
the spherical coordinate system. Under these conditions, the characteristics 
of droplet combustion and strand combustion approach one another and the 
results of Refs. 1 and 2 can be applied directly to the droplet combustion problem. 
Therefore, it is only necessary to convert the earlier zero-order results to the 
present notation. 


In Refs. 1 and 2, the expression for the liquid surface temperature is equi- 
valent to Eq. (17) and no change is required. Under the present assumptions, 
particularly unity Lewis number, the curvalinear coordinate system has no 
influence on the surface temperature. This is also true in cases where the Lewis 
number is not unity at the high activation energy limit where the reaction zone 
is collapsed to a flame sheet (11). 


In the case of strand combustion, 
terms of a parameter A^ as follows; 


A 

s 


*V* 


OOQ 


00 o 

in 


* 

(p V ) 
'^ooo ooo •' 


* 

C 

P 


the burning rate eigenvalue is given in 



* * * 
M R 

oo o 


(25) 


At the large droplet limit, with the combustion zone close to the surface, 
results 


. * 

m 

o 


*2 * * 
r p V 
s OOQ "O 


(26) 


Therefore, m^ is given by 

m = (A/A )^^2 (27) 

o s 

Knowledge of the droplet size and the chemical kinetic parameters provides A; 
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the quantity A is given as the eigenvalue of the zero order solution in Refs. 1 
and 2 , completing the determination of m^. 

2.4 Zero-Order Results and Discussion 

The physical property values used in the calculations were the same as those 
used in Refs. 1 and 2 and are representative of hydrazine decomposition. The 
reaction was taken to be second order, which correlates the bulk of the data in 
Refs. 1 and 2. The specific value of the various properties and kinetic parameters 
used in the calculations are summarized in Table 1. 

The first check of the zero order model involved comparing predicted and measured 
liquid surface temperatures as a function of pressure. As noted earlier, the 
theoretical prediction is independent of geometry under the assumption of unity 
Lewis number. Therefore, the comparison is identical to that made in Refs. 1 and 
2 between measured and predicted surface temperatures in the strand combustion 
case. These results have been replotted and are shown in Fig. 1. Th& agreement 
is good throughout the pressure range of the measurements (.51 - 19.8 atm). 

The next set of calculations involved determination of the steady state burning 
rates and temperature distributions. In addition to the general solution, the 
large and small A asymptotic cases (large droplet and small droplet limits) were also 
computed. In order to plot m^ as a function of A for the large droplet limit, 
values of A^ were obtained from the strand bombustion calculations. For low 
activation energies, A„ varies somewhat with pressure. However, for E* given in 
Table 1 A^ is relatively independent of pressure and the value 1.432 x 10' was used 
for all conditions (1,2). 

Some typical results of the calculations showing m as a function of A are 
given in Fig. 2. The conditions for these calculations involved adiabatic com- 
bustion at atmospheric pressure with an initial liquid temperature of 298K. The 
plot shows the small droplet limit (which is independent of A) ; the large droplet 
limit (where m^ is proportional to A); and the complete solution which provides the 
transition between these limiting regimes. In this case, the transition region falls 
in the range 10^ < A < 10^. 

Data on hydrazine droplet combustion was available from earlier work done in 
this laboratory. This involved droplet combustion in a combustion gas, under 
decomposition conditions for various droplet sizes at atmospheric pressure. 

The data Involved the ambient temperature range 1660 - 2530 K, and it was necessary 
to extrapolate the values to the 1345 K ambient temperature .of the present calcu- 
lations. The values used in the comparison in Fig. 2 are given in Table 2. The 
agreement between the predicted and measured burning rate is seen to be quite good, 
lending confidence to the theoretical model. Notably this correlation was obtained 
for a value of B* (in Table 1) which correlates the strand buring rate data at ele- 
vated pressures in Refs. 1 and 2. Thus it appears that the droplet measurements at 
atmospheric pressure are consistent with the second order reaction observed at 
elevated pressures for strands. 

Tenq)erature profiles from the complete solution are shown in Figs. 3 and 4 for 
the gas and liquid phase. As A Increases, the region of substantial temperature 
gradients is seen to approach the liquid surface. This manifests itself as a re- 
duction in curvature effects and approach to the large droplet limit where adiabatic 
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Table 1 

Properties used in the Theoretical Model 


Property 

Value 

a (N/m2) 

3.0078 X 10^^ 

B* (m^/Kg-s) 

8.385 X 10^^ 

Cl* (KJ/Kg-K) 

3.0961 

Cp* (KJ/Kg-K) 

3.0961 

E* (KJ/Kg-mol) 

111700 

L* (KJ/Kg) 

1715 

Lj* (KJ/Kg-mol) 

40794 

Le 

1 

M (Kg/Kg-mol) 

24 

n 

2 

q* (KJ/Kg) 

4937 

R* (KJ/Kg-mol-K) 

8.3144 

Ti* (K) 

298 


1345 

Y 

1.126 

6 

0 

Xj* (J/m-s-K) 

0.390 


0.176 

p^. (Kg/m^) 

1000 







100 
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Fig. 3 Steady state gas phase temperatures as a function of A at atmospheric pressure. 
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Fig. 4 Steady state liquid phase temperatures as a function of A at atmospheric 
pressure. 
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Table 2 

a 

Steady Droplet Burning Rate Data 


dg* (mm) 

it 6 

4iTm^ X 10 (Kg/s) 

1.3 

.70 

6.3 

1.4 

19.1 

135 


^Extrapolated to Tcoq = 1345K from data 
of Ref. 12 at atmospheric pressure. 


Table 3 

g 

Droplet Size Range of Transition Region 


P„* X lO"^ 

(N/m^) 

d * 

- d * ( y ) 


®min 

®max 

1 


25 - 

■ 7500 

10 


2.5 

- 750 

100 


.25 

- 75 


a 4 9 

Corresponding to 10 < A < 10 





-15- 


droplet combustion and strand combustion are essentially the same. 

Results at elevated pressures are similar to those at atmospheric pressure, 
with a range of 10^ < A < 10^ representative of the transition region where both 
curvature and reaction play a role In the steady droplet combustion reaction. 
Using properties In Table l,,.>the droplet size range for this range of A Is shown 
In Table 3. These ranges are somewhat broader than those suggested earlier 
(Ref. 7) based upon simple flame position arguments. In particular, the small 
droplet regime Is confined to smaller droplets at each pressure. Notably, the 
size range of technological Importance falls largely In the transition region. 


2.5 First-Order Solution 


2.5.1 First-Order Equations 


Utilizing the general zero-order results, the first order equations given 
In Ref. (7) reduce to the following. 

In the liquid phase, r < 1 


m^j^ = const. 

. o dT, m dT, , 2 m„- dT 

d ,2 1. _ _o 1 _ ^ fl o 

dr '• dr dr 6^ 1 dr 


(28) 

(29) 


In the gas phase, r > 1 


2 2 2 
T 3 -^ - iu)r T, = - Iwr T 
o dr 1 o 


(30) 


^ (r^ T , 
dr o dr 


d0 r lw0j 

) - m 


^ (r^ T 3 -^) - m 
dr o dr c 


dr 

T 


o 

dTi 


dr 

0, 

- T, 


— &) + 


= - (X^) ioj 


dT 


1 dr 
1 , 


E T 
T 


2 . n. 

r UJ T 


1 _ 


(31) 


‘^^o 2 ,Y- 1, , 2 

m, r (-' ) iw - r q w^ n 

1 dr Y o 


(32) 


The boundary conditions are; 
r = 0 


Tj^ = 0 


(33) 
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r = 1 


T = T, . 
fl 1 Is ls+ 


(34) 


dT dT dT 

^ '"o d7->s+ ^ - “o (35) 


d(e -T^) dT 

T ^ ) - T, x^) = m (6, - T,) + (1 - T 

o dr ' ^ Idr . o l 1 1 o 

s+ s+ 


- q) (36) 


L T 

0^ - = a q exp (-L^/T^) 


(37) 


Far from the droplet, 0^^ and must approach a constant value 


fiS 00 


h‘\‘ ■'i 


(38) 


The value of depends upon the case under consideration. For an unsteady gas 

phase, since 0, = 0 (t) at large r, the constant can be determined by considering 
the behavior of Eq.‘^(31) at large r. In this case, Kj=(y - 1/ y) which is the form 
for isentropic flow. 


In the case of a quasisteady gas phase, with transient effects still important . 
in the liquid phase, is the amplitude of the ambient temperature variation due 
to the varying energy content of the liquid at the surface of the droplet (re- 
sulting from transient energy storage in the liquid phase.) In this case, 
becomes a second eigenvalue in the solution to be determined along with m^. At the 
limit where the liquid phase also becomes quasisteady, the value of goes to zero, 
since in the absence of transient energy storage in the liquid phase, Eq. (11) 
must be satisfied. 

The equations for the first order are linear. Therefore, solutions can 
be combined to eliminate the need for iteration in order to determine mi and Ki. 
Since the limits of integration extend to infinity, asymptotic gas phase solutions 
are employed so that numerical computations can be confined to a finite region. 
Numerical difficulties are encountered in numerically integrating the liquid phase 
solution to r = 0, therefore an asymptotic analysis is employed in this region as 
well. 


2.5.2 First-Order Liquid Phase 

Similar to the zero-order case, the liquid phase solution is completed by 
applying the boundary condition T^^ = Tj^^ at r = 1, although T^^^ is not generally 
known until the gas phase solution is complete. This procedure allows the liquid 
phase solution to be decoupled' from the gas phase. 


In order to simplify subsequent discussion, Eq. (29) may be rewritten as: 


L {T^} 



H 


3 


(39) 



where 


-^"2 

dr 


is a linear operator, and 


0 

U =: _ • -Q 

r „ , 2 

PfS^r 

H = - ^ 

2 (S^ 


m 


m 


(Pffifr ) 


H = (T - T ) 5-2 exp (1 

3 os 1 r _2^'^ Pf°f 


7» 


(40) 


(41) 

(42) 

(43) 


The boundary conditions on Eq. (39) are 

r = 0, T = 0 ; r = 1 T = T, (44) 

1 * , 1 Is 

The quantity is an eigenvalue to be determined from the complete gas and 
liquid phase solutions. Exploiting the linearity of Eq. (39) let 

Ti=T^+mfiTi* (45) 

Equation (45) can be divided into homogeneous and particular parts 

■^1 ' "iH 

^1* - ’^IH 

^ * 

where K.J, and K.J, are constants to be determined from the boundary conditions. The 
differential equations then become 


L =0 ; L {Tj^p*} - H 3 (47) 

The values of K._ and iC, are selected to match the boundary condition at r = 1, 
to yield ^ 



-18- 


In order to match the liquid and gas phase solutions , the derivative of 
at the liquid surface is required. Combining Eqs . (45), (46) and (48), this 
quantity is 


di 


s- 


H, T, + m,, 

4 Is 5 fl 


(49) 


where 


dT 


in 


IHs- 


dr 


) . 


s- 


(50) 


«5 = 




) (-^) ] 


(51) 


s- 


The behavior of Eq. (47) at small r yields the following asymptotic forms 
for and T^^p* r ^ 0 



= exp {- m^/p^6^r} 
5^ «os - •'l> <1 



exp 


{ P- ■ 

Vf6f 



(52) 

(53) 


Equations (52) and (53) formally satisfy the boundary condition 0, as 

r 0. The arbitrary constants in homogeneous equation has been ^sorbed into 
^ and K^*. 

At the limit where = 0, the asymptotic forms of Eqs. (52) and (53) are 

valid throughout the range 0 < r < 1. At this limit, the solution becomes 


ui/6^ = 0 , 


m. 


[T 

^^Is ^ P.6 


rf 


'■^os - ■'iXi - 7>1 - 7» 


f'^f 


(54) 


The solution procedure involves employing Eqs. (52) and (53) to generate 
starting values for the numerical integration of Eqs. (47), at some small but non- 
zero value of r. The equations are integrated to the liquid surface, allowing the 
evaluation of and from Eqs. (50) and (51). Substitution of these quantities 
into Eq. (49) provides the needed expression in the gas phase boundary conditions. 
Once T^ and m., are determined from the gas phase solution, Eqs. (45), (46) and 
(47) allow complete determination of T^^ in the liquid phase. 

2.5.3 First-Order Quasisteady Gas Phase 

A number of simplifications are available for the analysis when the gas phase 
can be assumed to be quasisteady. This limit arises since frequency effects appear 
as 0 ) in the gas phase and ui/6^ in the liquid phase. At atmospheric pressure, for 
hydrazine, 6.” has the value 1715 and under these conditions transient effects 
in the liquid become important at frequencies substantially below the frequencies 
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where gas phase transients are encountered. Since 6^ is proportional to the mean 
pressure, the gap between the two characteristic frequencies progressively becomes 
smaller as the pressure increases. 


At the quasisteady gas phase limit, m^^ = mfi, a constant. For adiabatic 
combustion, consideration of the overall energy balance and integration of Eq. (31) 
shows that 9i = a constant as well. Equation (32) can then be simplified to 
yield 


where 


and 


L* {Tj^} = + A^ m^ + A^ 


^ ^ *1 I? ^^2 

dr 


, dX T m 

^1 ^ dr ^ r^ “ 
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r^T 
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X.2 rX dr X ^X “(1-X) ^2^ 

o dr o o o ' o' X 


- q n w 

^3 (1 - X^) 

dX 


A 
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r X 


o 

dr 


. - q w n 

A^ = ^ o 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 


Xhe boundary conditions at r ® 1 become 

^^1 

+ D„ X, - D, m, 
1 ci.r i 


where 


s+ 


'2 "Is *"3 '“1 


dX, 


^1 ^2 ^1 = ^3 ^1 ^ ^ 

s+ 


Vl = ""l ^2 


D, = X 
1 os 

dX 

■>2 ■ - “o (1 - B) + 

s+ 


®3 = ^os + ^ + »5 


(62) 

(63) 

(64) 

(65) 

( 66 ) 


(67) 
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F = - T 
1 os 


dT 

"2 = <-f ^ ^ ^ 

s- 


F- = m 
3 o 


a q L- 


= 1 + 


2 


G 2 = a q exp (- Lf/Tjjg) 
The outer boundary condition Is r -> <», 
Again exploiting linearity, define 


Ti=Ti+m^Ti +K^T^ 


( 68 ) 

(69) 

(70) 

(71) 

(72) 


(73) 


Separating the functions in Eq. (73) into homogeneous and particular parts 

/\ ^ ^ 


T* = + K* 


(74) 


^1 = "iP + ^ '^IH 


where 


r {T^^} = 0 


L' 




L' 


{T * } = 
^ IP 


L' 



(75) 


^ A 

The parameters K, K and K are constants which can be selected to satisfy any of 
the boundary conditions. Since the outer boundary condition can be handled 
by the asymptotic solution of Eqs. (75), Eq. (64) was chosen to determine the 
values of the K's. This yields 



lPs+ 


)/T 


lHs+ 


(76) 
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K = 


T * /.p 

lPs+' lHs+ 


(77) 
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The remaining two boundary conditions at r = 1, Eqs. (62) and (63), can be 
employed to determine the values of m^^ and . The result is 


^1 ^“^ll "^23 “ “^21 '^13^^ ^•^11 “^22 " -^12 “^21^ 


(79) 


where 


“l “ ^*^13 “^22 *^12 “^23^^ ^*^11 “^22 


dT 

J, T = D, + D„ T, * - D_ 

11 1 ' dr 8+ 2 ls+ 3 

dTj 

h2 = <dr 

dT 

“^13 “ " ^1 ^dF'^ + ~ °2 ^Is 
s+ 

dT/ 

•'21 = ^ (— ) ^ '2 ’'is* - ^ 

8 + 

•^22 ” ^1 ^dr ^ + F» T - F- 

s+ 2 ls+ 3 


Jl2 J21^ 


(80) 


(81) 


•^23 = ^ (dT) , - ^2 ^18+ 
s+ 


The derivatives in these expressions are obtained by differentiating Eqs. (74) 
noting that the K's are constants. 

In order to complete the solution, the asymptotic forms of Eqs. (75) 
at large r must be solved. Proceeding in the same manner as the zero-order 
case, the result is 


'^IH = ^IP = = 1 - exp (-i^/r) 

T^p = exp (-m^/r) 


(82) 


Equations (82) automatically satisfy the outer boundary condition of the gas phase 
solution. 

The solution is obtained by calculating D, , D2, D3, F^, F2, F^, and G2 
from the zero-order solution and the solution of the first-order liquid phase. 

The parameters are also available from the zero order solution. Equations 

(82) are used to generate starting values for the differential equations given by 
Eqs. (75), at some large but finite value of r. These^equations_are then inte- 
grated to the surface of the droplet. The parameters K,'K* and K can then be 
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determlned from Eqs. (76) - (78). With these parameters in hand, T, , T. , 

T- and the derivatives of these quantities can be determined from Eqs. (75) 
and the derivatives of these equations. Determining the from Eqs. (81) then 
allows K, and mj^ to be found from Eqs. 
evaluated as a function of r from Eq. 
from Eq. (45) for the liquid phase. 


(79) and (80). If desired, 

(73) for the gas phase^yielding T^^) and 


can then be 


2.5.4 First-Order Unsteady Gas Phase 

The methodology for the first-order unsteady gas phase solution is similar 
to that of the quasisteady gas phase. In this case m, and 0]^ are no longer constants 
and additional equations must be solved. However, the system of equations is 
linear so that solutions can be combined to provide a straightforward method without 
iterations. Asymptotic analysis provides starting conditions for the integration 
at a large but finite value of r, with integration proceeding to the liquid surface. 

A full description of this method will be deferred until the end of the next 
report period. 

As noted earlier, the strand results can be used directly at the large droplet 
limit. In order to consider response, the open-loop droplet response function 
is employed (8). In the present notation, this quantity is 

= Re {m. /m } (83) 

r 1 o 

where Re denotes the real part, i.e. that portion of the burning rate oscillation 
that is in phase with the pressure oscillation. For instability, must be 
positive and of order unity (the exact value depends upon the degree of damping 
present) at a point in the combustion chamber when the pressure is varying. 

The conversion of r^ and r^^, the strand burning rate and its perturbation, is 
straightforward at the large droplet limit, namely 

“l 

Re {^} = Re {r. } (84) 

m X 

o 

since r =1. The dimensionless frequency for the strand combustion case, 
is related to the present dimensionless frequency as follows: 


u_ = (A /A)w (85) 

S o 

As noted earlier, for an activation energy suitable for hydrazine, A = 1.432 x 
10^ at mean pressures in the range 1-100 atm. The above transformations are also 
valid in the quasisteady gas phase limit. 


2.6 First-Order Results and Discussion 

Calculations are still in progress for the quasisteady and unsteady gas phase 
response analyses.. Therefore, a discussion of the major results will be deferred 
until the next report period. In order to indicate the scope of the results, however, 
some initial findings of the analysis will be discussed here. As in the case of the 
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zero order analysis, properties and kinetic parameters suitable for hydrazine 
decomposition are employed for the numerical results (Table 1). 

Figure 5 illustrates results at the small droplet or evaporative limit, 
where reaction effects may be ignored. The open-loop droplet response function.^ 
and the magnitude of the surface temperature oscillation, Tj^g, are plotted as a \ 
f unction of pressure. The frequency for this plot is zero so that both the gas 
and liquid phases are quasisteady (The present analysis is not formally valid 
at 0 ) = 0, since the mean droplet size varies during an oscillation at this condition, 
therefore, the plot should be taken as representative of conditions where it is 
still valid to assume that both gas and liquid phases are quasisteady but with 
sufficiently high frequency so that changes in mean size can be ignored during a 
period of oscillation) . 


At this condition, the open-loop response function is negative, with a relatively 
small magnitude. This is in contrast to results for bipropellant combustion where 
the response function generally approaches zero or a positive limit for small u. 

With increasing pressure, the open loop response function becomes more negative. 


The surface temperature oscillation is in phase with the pressure oscillation, 
and its magnitude Increases as the pressure increases. This behavior results pri- 
marily from the form of the vapor pressure curve. At the quasisteady limit 


T 


Is ^ 



therefore since Tgg increases with increasing mean pressure, Tj^^ increases accordingly. 


At the same limit, using the result for T^^^ 

3 


- T 


m^ 7H 


os 


H * ■'os - '1 


This approximate formula indicates negative response at the quasisteady limit. With 
increasing pressure, T^g increases, providing a larger negative response. 


Physically, the negative response at the totally quasisteady limit is caused 
by the fact that the gasification rate decreases at the surface tenq)erature Increases 
for a constant ambient temperature. At the quasisteady limit, the ambient temper- 
ature (corresponding to the adiabatic flame temperature) is independent of temper- 
ature-satisfying this condition. During the rising portion of a pressure oscillation, 
the surface temperature increases in step to keep the mass fraction at the liquid 
surface approximately constant. This causes reduced rates of gasification at 
pressures above the mean and increased rates at pressures below the mean — providing 
a negative response. . 


- ' As the frequency is increased, within the quasisteady gas phase limit, the 
ambient temperature also begins to oscillate with finite magnitude. This is due to 
the thermal wave in the liquid phase varying the total energy content of the flow as 
it reaches the liquid surface. This effect compensates for the higher liquid tem- 
peratures at. pressures above the mean, tending to make the response positive. A 
second factor arises near the characteristic frequency of the liquid phase thermal 
wave. As the surface temperature begins to lag the pressure oscillation, theliquid 
temperature gradient at the surface becomes negative and supplements the energy flow 
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i-loop droplet response and surface temperature perturbation as a function of pressure 
he small droplet limit. 
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to the surface — increasing the gasification rate at pressures above the mean. 

This leads to a peaking of the response at this frequency level, but while the 
response is positive it is not very large (order 0.5). These results will be 
documented more fully in the next report. Based on present findings it appears 
that droplets at the evaporative limit do not contribute substantieilly to the 
response of the combustion process, due to their relatively low • response . 

Moving from the evaporative limit to the other end of the scale, the large 
droplet limit, much greater response is observed for monopropellant droplets. 

Figures 6, 7 and 8 illustrate the open- loop droplet response as a function of 
frequency for large droplets at mean pressures of 1, 10 and 100 atm. In addition 
to the quasisteady gas phase results, findings allowing for a completely unsteady 
gas phase are also shown. 

At the large droplet limit, the response is essentially unity at low frequen- 
cies for a second order reaction. With increasing frequency a peak is reached at 
frequencies near the characteristic frequency of the liquid phase. This arises 
largely due to an increase of flame temperature when the pressure is above the mean, 
resulting from optimum energy storage in the liquid phase. As might be expected, 
the magnitude of this peak depends upon the activation energy of the reaction, 
increasing with increasing activation energy. Beyond the peak, the quasisteady gas 
phase solution gives a continuously declining response, similar to the case of the 
evaporative limit. At higher frequencies, transient gas phase phenomena become 
Important, leading to a second peak in the open-loop droplet response function. 

Comparing Figs. 6, 7 and 8, it is evident that the magnitude of the first 
response peak increases with increasing pressure. Therefore, large droplets at' 
high pressures contribute most to the response, while small droplets contribute 
very little at all pressures. Response calculations in progress in the inter- 
mediate regime will complete information on the transition between these two limiting 
cases. 


2.7 Low-Frequency Limit and System Analysis 

As noted earlier, the previous response analysis is not valid at very low 
frequencies where the size of a droplet varies appreciably during a period of 
pressure oscillation. In order to obtain response in this region lifetime cal- 
culations must be undertaken, similar to those of Refs. 9 and 10, under conditions 
where an oscillatory pressure variation is Imposed upon the combustion process. 

Evaluation of the response of a spray at higher frequencies Involves a 
weighted combination of the response of all droplet sizes present in the combustion 
field. This calculation also requires computation of droplet life histories in 
order to determine the size distribution within the spray. 

Work on the formulation £ind analysis of these two related problems will begin 
in the next report period. The zero-order analysis for droplet response will be 
eiq)loyed to estimate the gasification rate of the droplets. This will not require 
any recalculation of zero-order results, since the existing findings (for example 
Fig. 2) can be correlated rather simply in terms of droplet size at any given pressure. 
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Fig. 6 Open-loop droplet response as a function of frequency at the large droplet limit 
and a pressure of 1 atm. 




Fi£. 7 Opeii-loop droplet respoase «s • function of frequency et the large droplet Unit 
and a pressure of 10 aca. 
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Fig. 8 Open-loop droplet response as a function of frequency at the large droplet limit 
and a pressure of 100 atm. 
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2 . 8 Summary 

An analysis has been developed to compute the combustion response of mono- 
propellant droplets (the open-loop droplet combustion response). Under conditions 
where the droplet size does not change appreciably during a period of oscillation, 
the method is developed as a perturbation analysis with the response given by 
the first-order solution. Asymptotic analysis is employed to treat boxmdary 
conditions at infinity and zero. 

Calculations were completed for the zero-order steady state at pressures of 
1-100 atm, and droplet sizes extending from the evaporative limit (where reaction 
effects are small) to the highly reactive limit where the combustion zone is 
close to the droplet. The temperature and burning rate predictions of this analysis 
were compared with earlier measurements of droplet combustion in an environment 
typical of an actual combustion chamber. The agreement between theory and experiment 
was excellent, indicating that the model gives a good representation of the com- 
bustion process. This provides some encouragement to continue with the response 
analysis. 

Response calculations are in progress. Work at the large droplet limit is 
complete for both a quasisteady and unsteady gas phase. Work at the evaporative 
limit and in the transition region for a quasisteady gas phase is nearing completion. 
Results for an unsteady gas phase under these conditions will be completed in the 
next report period. 

Analysis of the low frequency limit where the droplet size varies appreciably 
during a pressure oscillation, and the combination of individual droplet response 
to obtain the response of a spray are related problems. These calculations 
involve application of computations that are currently in progress and will be 
undertaken in the next report period. 


3. High Pressure Bipropellant Combustion 
3.1 Introduction 


Earlier work on the comparison between predicted and measured droplet temper- 
atures and gasification rates at elevated pressures in flowing combustion gases has 
been concluded. Results indicated satisfactory modeling at low pressures using a 
conventional model, and satisfactory modeling near the critical point using a high 
pressure phase equilibria model. 

Current effort involves development of an intermediate pressure model to provide 
a continuous transition between these two limiting cases. A second area of effort 
involved generalizing the results in order to provide a simple correlation for the 
evaluation of high pressure effects for given combustion chamber conditions. 


3.2 Intermediate Pressure Model 


The intermediate pressure model has been coded for the computer, but problems 
have been encountered in debugging the routine. In particular, difficulties have 
arisen in obtaining proper as 5 miptotic behavior toward the high pressure limit. 

Work is in progress to resolve these difficulties and results should be in 
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hand at the end of the next report period. 

Wilson parameters and Henry's constants have been accumulated for the 
paraffinic hydrocarbons and some of the air gases (N 2 , O 2 , H 2 O, CO 2 ) . The 
literature will continue to be surveyed in order to expand the data base of the 
analysis . 


3.3. Generalized Analysis 

Work with the generalized analysis has been deferred until the difficulties 
with the intermediate pressure analysis are resolved. 


3 . 4 Summary 

The intermediate pressure model has been coded,- but problems were 
encountered in obtaining the proper asymptotic behavior at elevated pressures. 
Work is in progress to resolve these difficulties and results should be available 
at the end of the next report period. 


i 
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